Phase diagram and dynamics of the projected SO(5)-symmetric model of high- Tc 

superconductivity 
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We present numerical studies of a quantum "projected" SO(5) model which aims at a unifying 
description of antiferromagnetism and superconductivity in the high-Tc cuprates, while properly 
taking into account the Mott insulating gap. Our numerical results, obtained by the Quantum 
Monte Carlo technique of Stochastic Series Expansion, show that this model can give a realistic 
description of the global phase diagram of the high-Tc superconductors and accounts for many of 
their physical properties. Moreover, we address the question of dynamic restoring of the SO(5) 
symmetry at the critical point. 
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SO (5) theory of superconductivity provides a unifying 
approach to describe both antiferromagnetism (AF) and 
superconductivity (SC) in the high-Tc cuprates 0. Mi- 
croscopic S0(5)-symmetric models [|[|^ successfully de- 
scribe many features of high-Tc cuprate physics and 
make a number of experimental predictions How- 
ever, the cuprates' Mott insulating behavior at half-filling 
severely challenges the validity of SO (5) theory. SO (5) 
symmetry requires collective charge-two excitations to 
have the same (vanishing) mass as collective spin-wave 
excitations and not the large Coulomb energy gap U of 
several eV observed in a Mott insulator at half-filling. 
To overcome this problem a low-energy effective model ~ 
the so-called "projected" SO (5) or pS0(5) model - has 
been introduced in which a Gutzwiller constraint 
of no-double-occupancy is implemented exactly. This is 
done by projecting out the particle pair excitations (the 
so-called 7r+ Goldstone mode) and by retaining only the 
massless magnon and hole-pair modes. The central hy- 
pothesis of pS0(5) theory is that this projected model 
is accurate in describing both the static and dynamic 
properties of the high-Tc superconductors (HTSC). For 
the fermionic excitations, a corresponding pS0(5) model 
has been shown to provide a natural explanation of the 
interrelation between SC and AF gaps [||. 

In this letter, we compare the properties of the pS0(5) 
model in two dimensions with a variety of salient features 
of the HTSC such as the global phase diagram and the 
neutron scattering resonance. The pS0(5) model turns 
out to give a good semiquantitative description of the 
SC phase, the behavior of the chemical potential /i as a 
function of doping and of the magnetic resonance peak. 
We identify a jump in the hole pair density S in the S{fi) 
diagram which is very similar to the available experi- 
mental data on La2-xSrxCu04 The jump signals 
coexisting hole-rich (metallic) and hole-free regions, the 
latter of which becomes antiferromagnetic at T = 0. We 
study the nature of the coexistence region, identify a tri- 



critical point where the coexistence line merges into a 
Kosterlitz-Thouless critical line, and discuss to what ex- 
tent the S0(5) symmetry can be restored at this point. 

All our results were obtained using a numerically exact 
Quantum Monte Carlo technique, the Stochastic Series 
Expansion (SSE) ]lO|-p^. SSE has three important fea- 
tures discerning it from other QMC techniques (see p^ ] 
for an overview). First, SSE is an exact method with- 
out any systematic error (such as the Trotter discretiza- 
tion error). Second, a non-local loop- update mechanism 
can easily be implemented within SSE |11|. This ap- 
proach, which is similar to the loop-algorithm |Q, al- 
lows for much lower autocorrelations times on large sys- 
tems or at low temperatures than any purely local update 
could reach . Such a favorable scaling behavior, which 
makes it possible to simulate systems with several thou- 
sand lattice sites, is neccessary for tracing phase tran- 
sitions of the Kosterlitz-Thouless type, as we encounter 
here. Third, we have developed a very efficient way to 
measure arbitrary Green's functions within SSE fl^ ]. 

We start with a more detailed presentation of the 
pS0(5) model. In Ref. Q it has been shown that the 
low-energy SO (5) excitations on the rung of a ladder can 
be cast into a picture of 5 hard-core bond bosons: three 
triplet states (<J,^2 3 4) ^'^d particle and hole-pair states 
(i^ and tj^). As an effective coarse-grained model, this 
description can be extended to a two-dimensional sys- 
tem, whereby the excitations are now defined on a 2 x 2 
plaquette [f 15|. The projection is implemented by re- 
stricting within the Hilbert space with tp{x) = 0. The 
projected SO (5) Hamiltonian takes the form 

H = A,Y^ tU^)t^{x) + {Ac- 2fi)^tl{x)t,,(x) (1) 

X .a = 2,3,4 X 

- Js ^ naix)na{x') ~ Jc^{tl{x)t,^{x') + h.C.) 
<xx' > ,a=2:3A <xx'> 

where Ua — {ta + iI)/\/2 are the three components 
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of the Neel order parameter. As and Ac (which mod- 
els the Mott-Hubbard gap U) describe the gap energy 
of the magnon and the hole-pair states at fJ, — 0, and 
the Js and Jc terms stand for the hopping and for cre- 
ation/destruction processes of these states. 

Let us now turn to the numerical results. We choose 
Js — Jc/'^, corresponding to the "S0(5)-symmetric" point 
determined in a mean- field analysis of the pS0(5) model 
0. J:=Js will be our unit of energy. Since the magnon 
density rapidly decreases as soon as one enters the su- 
perconducting region (see Fig. |l|), the precise value of As 
is not crucial for the dynamics of the model or the gen- 
eral structure of the phase diagram in this phase. This 
holds as long as As ^ 4 J, as for larger A there is no or- 
dered phase in mean field. We, thus, choose As = J, and 
shift the chemical potential so that Ac = As. The mean- 
field result predicts a phase transition from AF to 
SC accompanied by a jump in the hole-pair and magnon 
densities at /ic = 0, the state at fi — fic being a coherently 
mixed AF+SC phase without phase separation. How- 
ever, Gaussian fluctuations change the picture and pre- 
dict a first-order transition. Here, we want to study this 
region with an appropriate, i.e. exact numerical evalua- 
tion. In fact, we expect the picture to change appreciably, 
since no long-range order is allowed in two dimensions. 
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FIG. 1. Hole concentration 5 ^ p/2 = \{t\t,^) (filled 
squares) and magnon density ^ ^{tata) (circles) as a func- 
tion of the chemical potential /i at T/J = 0.03. The plot- 
ted points result from a finite-size scaling with lattice sizes 
1/ = 10 X 10, 14 X 14 and 20 x 20. The small inlay shows a 
detailed view to the fi region in which the hole-pair density 
jumps to a finite value. The additional solid lines with error 
bars are T = data obtained from a simulataneous scaling 
of /3 oo and V -+ oo (with lattice sizes of F = 8 x 8, 10 x 10, 
12x12, 16x16, 20x20 and /3 = 4.8, 7.5, 10.8, 19.6, and 30). 

In Fig. 0we plot the mean hole-pair and magnon densi- 
ties as a function of the chemical potential for T/J = 0.03 
and their T— >0 extrapolations. The jump can be clearly 
seen at fic = —0.175, shifted with respect to the mean- 
field value due to the stronger fluctuations of hole pairs 
as seen in the Gaussian contributions j^]. The size of 



the jump, Sc « 0.11, as well as as its slope for fj, > ^c, 
9(5/9/i « 0.13, are very close to the mean-field values of 
1/8 [|l6|. Notice that, while the hole-pair density rapidly 
vanishes below /i^, a nonvanishing magnon density per- 
sists well beyond /ic even at T = 0, due to the fact that 
the magnon number is not conserved. 

A jump in the density as a function of /i has also been 
measured in La2_xSrxCu04 ||^. A comparison with the 
experimental data of Ref. |^ is displayed in Fig. ^ [ p^ . 
As one can see, our data reproduce the experimental re- 
sults within error accuracy. The data are fitted by ad- 
justing the energy scale J to 220 meV. This value of J 
has the correct order of magnitude for a typical magnetic 
superexchange interaction in the cuprates p5[ . 
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FIG. 2. Chemical potential shift A/i as a function of hole 
concentration 5 (open circles) for La2-xSrxCu04 at T = 77 
K. The diamonds and the solid line show numerical data for 
the infinite lattice obtained by finite-size scaling with J — 220 
meV and T/J = 0.03, which corresponds to T~77 K. 
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FIG. 3. Hole densities of the coexisting phases on the first 
order transition line from (almost) zero to finite hole density 
at /i = /Xc as a function of temperature. 

The nature of the phase transition at fJ. = fic can be de- 
termined by studying histograms of the hole-pair distri- 
bution for fixed fi = ^c- While in an homogeneous phase 
the density is peaked about its mean value, at fi = fic 
we obtain two peaks indicating a first-order transition 
with a phase separation between (almost) hole-free re- 
gions and regions with high hole-pair density. From Fig. 
^ we see that the transition is first order for T <Tp — 
(0.20 ±0.01) J at /x = ^p = (-0.168 ±0.002) J. Above Tp, 
the histograms show strongly fluctuating hole-pair densi- 
ties suggesting the presence of critical behavior. Indeed, 
as we will show below, the Kosterlitz-Thouless transition 
line between the SC and the paramagnetic phase merges 
at the point P, which is thus a tricritical point. 
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In realistic electron systems, however, long-range 
Coulomb repulsion between the doubly-charged hole 
pairs should heavily disfavor phase separation. To study 
the effect of off-site Coulomb interaction we have added 
additional nearest-neighbor and next-nearest-neighbor 
Coulomb repulsions V and V = 0.67 V to the pS0(5) 
model. Indeed a relatively moderate Coulomb repulsion 
of y/ J « 0.2 is enough to completely destroy the phase 
separation. The interesting effect of Coulomb interac- 
tion in two dimensions is, thus, to push down the tricrit- 
ical point into a quantum-critical point at T — Q. This 
quantum-critical point now separates the SC from the AF 
phase and could be a possible candidate for the dynamic 
restoring of SO (5) symmetry. 

The next question concerns the nature of the phase 
above the tricritical point (or above the quantum critical 
point in the presence of Coulomb repulsion): do we find 
microscopic density modulations in the form of static or 
dynamic stripes or a phase which is homogeneous also 
on the microscopic level? Careful studies of the static 
density-density correlation (n/i(r)n/j(0)) and the dynam- 
ical response function D(k, w) = (n?i(k, w) n/i(k, 0)) show 
that there is no sign of stripe formation, neither without 
nor with the above (moderate) Coulomb repulsion. 

We have anticipated the existence of a SC phase for 
/i > /ic- In fact, at T > a true long-range order is pro- 
hibited by the Mermin- Wagner theorem. However, we 
can still have a KT phase of finite superfluid density ps 
at finite temperature, which is identified by a power-law 
decay of the SC correlation function 

Ch{r) = {t\{r) + t^(r)) (^(O) -|- ^^^(0)) . 

The transition separates long-range power-law {Ch{r) oc 
r~") from rapid exponential decay (C/i(r) oc e"^*"). A re- 
liable and accurate distinction between these two decay 
behaviors requires a finite-size scaling with large system 
sizes, as well as an efficient QMC estimator for the Green 
functions appearing in the correlation function. With its 
non-local update scheme and with our new estimators 
for arbitrary Green functions SSE provides both. An al- 
ternative method for detecting a KT transition exploits 
the fact that the superfiuid density jumps from zero to 
a finite value at the KT temperature Tkt [Q- Since 
within SSE the superfluid density can be measured quite 
easily by counting winding numbers , this criterion is 
numerically even more favorable than the arduous direct 
determination of decay coefficients. 

Fig. 1^ plots the phase diagram obtained by apply- 
ing both criteria independently. The figure shows that 
the projected S0(5) model indeed has a KT phase with 
quasi long range order whose form in fi-T space looks 
like the one of the high-Tc cuprates. Both criteria re- 
sult in exactly the same clearly pronounced phase sepa- 
ration line. It is well known that a similar transition can- 
not occur for antiferromagnets , and that the finite-T 



AF correlation length ^ is always finite and behaves like 
^ oc Q'^^Ps/knT ^ being the spin stiffness. This fact is 
confirmed by our numerical results. 

In the preceeding paragraphs we have presented some 
general properties of the pS0(5) model. The question 
"how S0(5)-symmetric is the pS0(5) model", i.e. the 
question whether there exists a point on which the full 
S0(5)-symmetry is dynamically restored, still remains 
open. As one can see from Eq. (|l|), the excitation en- 
ergy for hole pairs can be compensated by /z, in order 
to have equal energies for local spin and hole-pair ex- 
citations. Due to this partial compensation the mean- 
field ground state of this model recovers exact S0(5) in- 
variance at Jc = '2, Jg and = Ac 0. However, since 
the Casimir operator of the SO (5) group does not com- 
mute with the Hamiltonian, this invariance is not exact, 
and a symmetry breaking effect can already be seen at 
the Gaussian level |^]. For a classical three-dimensional 
S0(5)-symmetric model, numerical simulations indicate 
that the symmetry is asymptotically restored at a bicrit- 
ical point provided the symmetry-breaking terms have 
the appropriate sign ||5|,pO|. This is in contrast with the 
prediction from the e-expansion pH , which would sug- 
gest a fluctuation-induced first-order transition. On the 
other hand, other arguments seem to indicate that the 
decoupled fixed point is slightly attractive |22|. This 



discrepancy clearly indicates that strong- coupling effects 
play an important role, calling for an appropriate, i.e. 
numerically exact treatment, as it is provided by SSE. 

One neccessary condition for an S0(5)-symmetric 
point is that the formation energies of hole-pair bosons 
and of magnons are identical. This condition is fulfilled 
along the line from S to the tricritical point P in Fig. 
^. Another neccessary condition is that hole pairs and 
magnons behave in the same way at long distances. This 
condition is fulfilled on the dashed line in Fig. ^, where 
the AF and SC correlation lengths ^ become equal. In- 
terestingly, these two conditions meet (within error bar 
accuracy) at the tricritical point P. Of course, here the 
correlation length is still finite; however, we find rela- 
tively large f values of order 10 to 15 in the immedi- 
ate vicinity of point P, demonstrating the importance of 
SO (5) critical fluctuations in this region. 
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FIG. 4. Phase diagram of the pSO(5) modeh The squares 
between S and the tricritical point P trace the first-order line 
of phase separation. The solid line from P to the right edge 
of the plot traces the Kosterlitz-Thouless transition between 
the "SC" and the normal state. The dashed line separating 
Nt (^triplet dominated region) and Nh (=hole pair domi- 
nated region) describes the line of equal AF and SC correla- 
tion lengths. The small inlay shows the same phase diagram 
on a larger /i scale, covering the whole "SC" KT phase. 

One of the main features of S0(5) theory is that it pro- 
vides an elegant explanation for the neutron resonance 
peak observed in some high-Tc cuprates at fc = (tt, tt) [Q . 
Experiments show that the resonance energy Wros is an 
increasing function of Tc, i.e. Wres increases as a func- 
tion of doping in the undcrdoped and decreases in the 
overdoped region |23|. Fig. ^ shows the corresponding 
spin correlation spectrum obtained from the projected 
SO (5) model. The "AF" and underdoped SC regions 
are described correctly: spin-wave excitations are mass- 
less Goldstone modes in the AF phase at fi < fic (and 
T = 0), and become massive when entering into the SC 
phase. Wres increases monotonically up to optimal doping 
/^opt ~ 1. In the overdoped range, however, Wres is further 
increasing in contrast to what happens in the cuprates. 
This is not surprising given the fact that SO (5) theory 
was developed to model the interplay of antiferromag- 
netism and superconductivity in the vicinity of the AF- 
SC transition, i.e. in the underdoped range. The reso- 
nance peak continuously looses weight when increasing fi, 
which is consistent with experimental observations l23] . 
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FIG. 5. Dispersion of the (7r,7r) peak of spin correlation as 
a function of the chemical potential The numbers in paren- 
theses indicate the peak weights, i.e. the area under the peak. 
(20x20 lattice at temperature T/J^O.l). 

A comparison of the critical temperature Tc obtained 
from Fig. ^ and Wi-os at optimal doping yields the ratio 
Tc/ujycs,opt = 0.23. This is again in accordance with the 
corresonding ratio for YBa2Cu306+x, for which the ex- 
perimentally determined values Tc = 93 K (thus ksTc = 
8.02 meV) and Wros,opt = 41 meV yield Tc/wios.opt = 0.20. 

In summary, we have shown that the projected SO (5) 
model - which is the simplest bosonic model contain- 
ing two key ingredients for the HTSC: the Mott gap and 
the vicinity of AF and SC phases - gives a semiquanti- 



tative description of many properties of the HTSC in a 
consistent way. In particular, we have identified a KT- 
SC phase whose bounds look similar to the the ones of 
real cuprate materials. Also the doping dependence of 
the chemical potential as well as that of the neutron res- 
onance peak in the underdoped regime are reproduced 
correctly. Although the second issue discussed in this 
letter, the dynamic restoring of SO (5) symmetry in the 
vicinity of a critical point, cannot be conclusively an- 
swered in D = 2 dimensions, we have found a critical 
point scenario which makes a partial dynamic restoring 
of S0(5) symmetry at a T ^ bicritical point in _D = 3 
a good possibility. This will be discussed in a future 
publication. The authors acknowledge financial support 
from BMBF (05SB8WWA1), DFG (HA 1537/16-1,2 and 
HA 1537/17-1), KONWIHR OOPCV and the Swiss Na- 
tional Science Foundation. The calculations were carried 
out on a Cray T3E at the high-performance computing 
center HLRZ (Jiihch). 
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